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Collimation and Radiative Deceleration of Jets in TeV AGNs 
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We consider some implications of the rapid X-ray and TeV variability observed in M87 and the TeV blazars. We outline 
a model for jet focusing and demonstrate that modest radiative cooling can lead to recoUimation of a relativistic jet in 
a nozzle having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances 
from the central engine and may explain recent observations of the HST-1 knot in M87. Possible applications of this 
model to TeV blazars are discussed. We also discuss a scenario for the very rapid TeV flares observed with HESS and 
MAGIC in some blazars, that accommodates the relatively small Doppler factors inferred from radio observations. 
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1. Introduction 

Recent observations of the VHE emission from com- 
pact relativistic systems raise new questions regard- 
ing the dynamics and dissipation of relativistic jets. 
Of particular interest is the class of TeV AGNs for 
which X-ray and TeV observations indicate: (i) ex- 
tremely rapid variability of the VHE emission, in two 
cases over timescale significantly shorter than the dy- 
namical time Tg/c of the putative, central black hole; 
(ii) substantial differences between the Doppler fac- 
tor inferred from TeV observations and that associ- 
ated with radio knots; (iii) rapid variations of the flux 
emitted from a region located far (at radii > lif'rg) 
from the central engine. These observations motivate 
reconsideration of the standard view, according to 
which dissipation of the bulk energy on sub-parsec 
scales is accomplished predominantly through forma- 
tion of internal shocks in colliding fluid shells. Below 
we review the observational motivation in greater de- 
tail and then go on to discuss the implications for the 
collimation and dissipation of the outflow. 

1.1. TeV blazars 

Large amplitude variations of the VHE 7-ray flux 
on time scales of minuets have been reported for 
Mrk 421, Mrk 501 and PKS 2155-304. The im- 
plied size of the region producing these 7-ray flares 
(as measured in the Lab frame) is limited to 
Ar ^ 10^'*r(5tvar,h/ (1 + ^) cm, but is naively expected 
to exceed r^. Here F and 5 are the bulk Lorentz 
factor and the corresponding Doppler factor of the 
emitting matter, respectively, tvar,h is the observed 
variability time in hours, and z is the redshift of the 
source. The extremely short variability reported re- 



cently for PKS 2155-304 [1] and Mrk 501 [2] is there- 
fore puzzling, as it implies that either, the black hole 
mass Mbh ^ lO''Af0, inconsistent with other esti- 
mates [3], or that the source imprinting the variabil- 
ity has a characteristic size considerably smaller than 
the black hole's horizon. The requirement that the 
pair production opacity at TeV energies is small im- 
plies high Doppler factors, (5 ~ 30 — 100 [4, 5], if the 
TeV emission originated from jet radii Tom ~ Ar. 
Such high values are consistent with those obtained 
from fits of the SED to a homogeneous SSC mode, 
but are in clear disagreement with the much lower 
values inferred from unification schemes [6, 7] and 
superluminal motions on parsec scales [8-10]. 

It has been argued that such high values of the 
Doppler factor may not be required if the 7-ray pro- 
duction zone is located far from the black hole, at 
radii rem >> Ar. This possibility is motivated by 
recent observations of M87, as described below. In 
that case the compactness of the TeV emission zone 
may be constrained by the variability of the IR flux 
observed simultaneously with the TeV flare, allowing 
low values of S in cases where the variability time 
of the IR emission is much longer than the dura- 
tion of the TeV flare. However, such a mechanism 
requires either jet power much larger than the lu- 
minosity of the TeV emission measured during the 
flare, LxeV 10^^ erg s^^, or focusing of the jet. 

1.2. M87 

The HST-1 knot in M87 is a stationary radio fea- 
ture associated with the sub-kpc scale jet. The knot 
is located at a projected distance of 60 pc (0.86") 
from the central engine, and is known to be a re- 
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gion of violent activity. Sub-features moving away 

from the main knot of the HST-1 complex at su- 
perluminal speeds have been detected recently [11]. 
In addition, rapid, large amplitude variations of the 
resolved X-ray emission from HST-1 have been re- 
ported, with doubling time tvar as short as 0.14 yrs. 
The observed variability limits the linear size of the 
X-ray source to Ar ^ TSty^T ~ 0.022r^ pc, which for 
reasonable estimates of the Doppler factor is three 
orders of magnitude smaller than the distance be- 
tween the HST-1 knot and the central black hole. 
Based on a claimed correlation between the X-ray 
and TeV emission it has been proposed that HST- 
1 may also be the region where the TeV emission is 
produced [11] (but see, e.g., Ref 12 for other explana- 
tions). As mentioned above, this motivated the con- 
sideration that the TeV emission zone in TeV blazars 
may also be located far from the black hole. 

It has been proposed that HST-1 reflects the 
location of a rccoUimation nozzle [13]. In this pic- 
ture the superluminal sub-knots that seem to be ex- 
pelled from the HST-1 complex can be associated 
with internal shocks produced by reflection of the 
recollimation shock at the nozzle. The rapid vari- 
ability set a limit on the cross-sectional radius of the 
channel at the location of HST-1 that depends on 
the fraction of jet power radiated as X-rays (and 
TeV emission, if indeed originating from the same 
location). Estimates of the jet power in M87 yield 
Lj > 10^4 erg s"^ [13, 14]. For the observed X-ray 
[11] and TeV [15] luminosities, LxeV ~ ~ 10''^ 
erg s^^, this implies a rather small radiative effi- 
ciency, 7] = Lx/Lj ~ 10 In order to account for 
the observed luminosity the scale of the fluctuations 
producing the X-ray flare must satisfy d ^ rj^^^a, 
where a denotes the cross-sectional radius of the 
jet at the dissipation region (see Fig. 1). The vari- 
ability time, on the other hand, constrain the size 
of these fluctuations to he d ^ FiJ^var. Combining 
the two constraints yields a ^ ry~^/^r^fvar- Associ- 
ating the apparent speed measured for the super- 
luminal sub- knots in the HST-1 complex with the 
Lorentz factor of the fluid in the vicinity of the noz- 
zle gives r ^ 4 [11]. Adopting F = 4 and viewing an- 
gle 6n = 30°, we estimate a ^ 1 pc for the reported 
X-ray variability. Assuming that the TeV emission 
is associated with HST-1 would require an even bet- 
ter focusing, a ^ 0.1 pc. At the radius of HST-1 this 
implies a/rnsTi ~ 10"'^. 




Fig. 1. Sketch of a rccoUimation nozzle. The shaded region 
represents an emitting blob of size 2d. 

It should be noted that the TeV luminosity ob- 
served typically in TeV blazars, LxeV ~ lO^''"^^ gj-g 
s~^ , is much larger than X-to-7 ray luminosity seen in 
M87, and is likely to be comparable to the jet power. 
If the rapid TeV variability in the TeV HBLs is to be 
explained by recollimation of a mildly relativistic jet 
(with fluid Lorentz factor of the order of the apparent 
superluminal speeds measured in HBLs), then this 
would imply far more stringent constraints on the 
nozzle. Whether the extreme TeV flares observed in 
VHE blazars can be accounted for by recollimation 
shocks at radii rem >> Ar remains to be explored. 
Recollimation shocks may be an important dissipa- 
tion channel also in other sources, e.g., GRBs [17]. 

Below we demonstrate that even modest radia- 
tive cooling of the shocked jet layer deflected by the 
external medium can lead to extremely good focus- 
ing of the channel. This has been studied earlier in 
the non-relativistic regime in the context of SS433, 
e.g., Ref. 16 

2. Focusing of a Radiative Flow by 
External Medium 

Bromberg and Levinson [17] (hereafter BL07) con- 
structed a class of semi-analytical models for the 
confinement and coUimation of a relativistic jet by 
the pressure and inertia of a surrounding medium. 
They considered both, confinement by kinetic pres- 
sure of a static corona, and confinement by the ram 
pressure of a supersonic wind emanating from a disk 
surrounding the inner source. In general, the colli- 
sion of the inner jet with the confining medium leads 
to the formation of a contact discontinuity across 
which the total pressure is continuous, and an oblique 
shock across which the streamlines of the colliding 
flow are deflected. In cases where confinement of the 
inner flow is accomplished through collision with a 
supersonic wind a second shock forms in the exte- 
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Fig. 2. Lorentz factor of the shocked and unshocked jet in the 
case of confinement by a static corona with a pressure profile 
Pext z~^, for different values of the cooling parameters. The 
ratio of the total luminosity radiated away by the shocked jet 
layer and total jet power is Lc/Lj = (no radiative losses) 
in the upper panel, Lc/Lj = 0.1 in the middle panel, and 
Lc/Lj = 0.27 in the lower panel. The injected flow in all 
panels consists of a cold, purely baryonic fluid with Lorentz 
factor r = 10 at the injection point (2 = 0). The shock surface 
is marked by the dashed line 

rior wind. The model outlined in BL07 computes the 
structure of the shocked layers of the deflected inner 
jet and the exterior wind in the latter case, assum- 
ing a steady, axisymmetric flow. In BL07 the focus 
was on the application to GRBs. Radiative losses 
have been ignored since the large optical depth of 
the shocked jet layer on scales of interest renders 
such losses negligibly small. In blazars, recollimation 
shocks are expected to form above the photosphere. 
If a non-negligible fraction of the energy dissipated 
behind the shock can be tapped for the acceleration 
of electrons to nonthermal energies then the cooling 
rate may be large enough to affect the structure of 
the recollimation nozzle. To study this effect we in- 
corporated synchrotron colling into our model. To 
be more specific, we assume that a fraction of the 
total energy flux is carried by magnetic fields, and 
that a fraction of the enthalpy behind the shock 
is injected as a power law distribution of electrons: 
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Fig. 3. The effect of jet composition on the coUimation pro- 
cess. In all cases shown Lc/Lj = 0.27. The ratio of total jet 
energy to rest mass energy is Lj/MpC^ = 1 in the lower panel 
(corresponding to the cold, baryonic jet exhibited in Fig. 2), 
Lj/MpC^ = 4 in the middle panel, and Lj /MpC^ = 10 in the 
upper panel. 

dUe/dte oc e^^. 

Examples are shown in figures 2 and 3. In all 
the examples shown the confining medium is static 
with a pressure profile p^xt oc z^^. In Fig. 2 the jet 
content is dominated by baryons; that is, at the in- 
jection point Lj ~ MpTc? (the energy density of the 
magnetic field in all examples shown here is taken 
to be at most a few percent of the total energy den- 
sity). The difference between the three cases shown 
is in the synchrotron cooling rate behind the oblique 
shock (which we control by the parameters and 
^s). The ratio of the total luminosity Lc radiated 
away by the shocked jet layer and total jet power Lj 
injected at z = in the upper, middle and lower pan- 
els is Lc/Lj — (no radiative losses), 0.1, and 0.27, 
respectively. As seen, in all cases the shock (indicated 
by the dashed line) converges to the axis at roughly 
the same distance from the injection point. The main 
effect of the radiative cooling, as seen in Fig. 2, is to 
increase the shock compression ratio, thereby reduc- 
ing the width of the shocked layer and, as a result, the 
cross-sectional radius of the jet at the recollimation 
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nozzle. Reflections of the converging shocked layer at 

the symmetry axis should lead to formation of inter- 
nal shocks in the vicinity of the recoUimation noz- 
zle, as seen in full 2D simulations. We propose that 
the superluminal motions of the sub-features in the 
HST-1 complex reported recently are associated with 
these internal shocks. A considerable fraction of the 
remaining jet power can dissipate via these internal 
shocks in a region much smaller than what would be 
expected in the case of a conical jet. This may lead 
to large amplitude variations over timescales much 
shorter than the jet radius. 

The effect of cooling is less dramatic in light jets 
in which the specific enthalpy of the unshocked jet 
is larger than unity, as demonstrated in Fig. 3. The 
reason is that in that case the bulk energy is domi- 
nated by pressure rather than rest mass energy and, 
therefore, significant fraction of the jet energy must 
be radiated away in order that the shocked jet layer 
will become under-pressured. The conditions under 
which eS'ective focusing occurs are under investiga- 
tion. 

3. Radiative Deceleration of Fluid 
Shells as the Origin of Rapid Flares 
in TeV Blazars 

An alternative scenario for the rapid TcV flares has 
been proposed recently [18], in which flares observed 
in sources like Mrk 421, Mrk 501 and PKS 2155- 
304 are produced by radiative deceleration of fluid 
shells expelled during violent ejection episodes. This 
model cannot naturally account for flares of dura- 
tions considerably shorter than Tg/c, as discussed 
above. The add hoc assumption made here is that 
fluid shells accelerate to a Lorentz factor To >> 1 at 
some radius ~ 10^ — lO^r^, at which dissipation 
of their bulk energy occurs. The dissipation may be 
accomplished through formation of internal shocks 
in a hydrodynamic jet or dissipation of magnetic en- 
ergy in a Poynting flux dominated jet [19, 20], and it 
is assumed that a fraction of the total proper jet 
energy density, ti'^ , is tapped for acceleration of elec- 
trons to a maximum energy jmaxTneC^ ■ The shocks 
may also result from a focusing of the outflow, as de- 
scribed in the previous section, which can in principle 
give rise to variability over timescale of the order of 
oF^^ < Tg/c, where a again is the cross-sectional ra- 
dius of the nozzle and r„ the Lorentz factor of the 
flow passing the nozzle. Whether the pressure in the 



ambient medium is sufficient for significant focusing 
of the jet and on what scales is yet an open issue. 
The fluid equation. 



^ rpOfJ. 



(1) 



where S*^ accounts for energy losses due to radiative 
friction, can be solved, assuming that the intensity 
of ambient radiation depends on radius as /« oc 
and that the proper density and average energy of 
the nonthermal electrons are independent of radius, 
to yield the asymptotic Lorentz factor [18]: 



I 



Too^roy^, (2) 

where Fq = r(r = r^) is the initial Lorentz factor. 
The stopping length I can be expressed in terms of 
the optical depth for 77 absorption of a 7-ray of en- 
ergy rrieC^e-y by a power law target photon fleld of 
the form Is{is) « e^"; eg, mm < < ts,max-, as: 

- = i^^) 9{e,), (3) 

where ^(e^) < 1 is some function of energy that de- 
pends on the spectral index a, and is given explic- 
itly in Ref 18, and x =< 7^ > /(< 7 > 7max) is a 
dimensionless factor that depends on the energy dis- 
tribution of nonthermal electrons. For a power law 
distribution, dng/d'y oc 7"' with q < 2, we have 
1 >X> 0.1. 

The main conclusion from Eq. (3) is that for a 
reasonably flat distribution of nonthermal electrons, 
q < 2, extension of the distribution to a maximum 
energy 7max at which the pair production optical 
depth, r(ro7max), is a few is already suflacient to 
cause appreciable deceleration of the front. It can be 
shown [18] that for the TeV blazars a background 
luminosity of Lg ~ 10^^ — 10*^ erg s~^, roughly 
the luminosity of LLAGN, would lead to a substan- 
tial deceleration of the front and still be transparent 
enough to allow the TeV 7-rays produced by Comp- 
ton scattering of the background photons to escape 
the system. The ambient radiation field is most likely 
associated with the nuclear continuum source. Prop- 
agation of the 7-ray flare from low-to- high energies, 
as reported recently for Mrk 501 [2], are naturally 
expected in this model, since the 7-spheric radius 
increases quite generally with increasing 7-ray en- 
ergy [21]. The bulk Lorentz factor of the jet dur- 
ing states of low activity may be appreciably smaller 
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than that of fronts expelled during violent ejection 
episodes. 
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